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' \ Abstract 

Using finite-size scafing methods we measure the thermal and mag- 
• netic exponents of the site percolation in four dimensions, obtaining 

a value for the anomalous dimension very different from the results 
found in the literature. We also obtain the leading corrections-to- 
''^ • scaling exponent and, with great accuracy, the critical density. 
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1 Introduction 



From the point of view of its definition, the simplest statistical system is 
perhaps the percolation. In the case of the site percolation, we fill the sites 
of a given lattice with probability p. Then we construct the clusters as sets 
of contiguous filled sites. 

The critical properties of the system can be described in terms of the clus- 
ters. For instance, at the critical percolation the mean cluster size diverges. 
We define the percolating cluster as the one that contains, in the thermody- 
namical limit, an infinite number of sites. The strength of this cluster (i.e. 
the probability of containing an arbitrary point) is the order parameter of 
the transition: it is zero for p < Pc, and finite for p > Pc 

Another interesting model is the bond percolation. In this case we fill the 
lattice bonds with a given probability and construct clusters analogously. It 
is believed that both models belong to the same Universality Class (share 
the critical exponents). 

It is possible to relate the percolation problem (in the bond version) with 
the g-states Potts model using the "Fortuin-Kasteleyn" representation of the 
latter. The bond percolation is obtained in the q —>■ 1 limit 

Moreover it is possible to write down a field theoretical description of the 
percolation. In general, the Potts model is described by means of a 0^ theory, 
where the coefficient of the cubic term is proportional to g — 2. For the Ising 
model {q = 2) this term vanishes, and the leading term is 0"^, recovering the 
standard field theory representation. For q ^ 2 we can write 



S= d'^x 
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-(V0i)(V</)i) + -ml(j)i(j)i + -^gQdijk(t)i(j)j(l)k 
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where the coefficients d^k depend on the model (Potts, percolation, Lee- Yang 
singularities, etc.), and n = q — 1 is the number of components of the field 
0i. Thus, the percolation is described by the action (|^) in the limit of zero 
components of the fields. 

Using the standard tools it is possible to obtain an e-expansion for this 
model (and in particular for the percolation). The power counting tells us 
that the upper critical dimension of the model is six and thereby the expan- 
sion parameter is e = 6 — Results up to three loops can be found in the 
literature 0. 

For large dimensions {d = 5, and, of course, 6) there is a good agreement 
between the results obtained from the e-expansion (resumed using Pade tech- 
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niques), the values from numerical simulations, and the results from high 
temperature expansions. 

In lower dimensions, the results disagree for the anomalous dimension, 77. 
The e-expansion predicts a clear negative value, while in the two dimensional 
case rj should be non-negative because the correlation function is decreasing 
with the distance. In fact, in this case, it has been conjectured 0] that 
r] = 5/24. 

In this paper we will show that the value of the four dimensional 77 expo- 
nent turns out to be large by a 30%, compared to the e-expansion. Thereby 
it remains as an open problem to understand why the convergence of the 
e-expansion for this model is so poor even for small values of e [Q. In order 
to calculate critical exponents we extend some recently developed accurate 
finite-size scaling techniques to site percolation. As a benchmark we re- 
port the two dimensional critical exponents (for which there are almost exact 
analytical estimates). 

A related model with the site percolation is the diluted Ising model 0. 
It is defined as a standard Ising model where the spins live only on filled 
(with probability p) sites. The field theoretical description of this model is a 
^''-theory with a random mass term. Using the replica trick it can be related 
with an O (N) symmetric 0^ theory with cubic anisotropy, in the limit of zero 
field components (i.e. 0) 0. 

The limit of zero temperature (large (3) of the diluted Ising model is the 
site percolation while when p ^ 1 it is the pure Ising model. A precise 
determination of the critical exponents of the d = 4 percolation is also a very 
useful first step to understand the phase diagram (/3,p) in the diluted Ising 
model. On the other hand, the site percolation is useful as a benchmark to 
develop and test different tools to apply to more complicated systems as the 
d = 4 diluted Ising model [0. 

Finally, we remark that we are specially interested in these four dimen- 
sional models in relation with the triviality issue (is there an interacting 
continuum limit in four dimensions?). In order to solve the triviality prob- 
lem is crucial to characterize all the possible fixed points in four dimensions. 
The site percolation has the unusual feature of having the critical dimension 
at d = 6, thus, it does not present the usual Mean Field exponents at = 4. 
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2 Numerical Methods 



We will work in a liypercubic lattice of linear size L with periodic boundary 
conditions. The Monte Carlo (MC) procedure for generating configurations 
in this model is straightforward: we fill each lattice site with probability p. 
The next step is to build the clusters, what is a deterministic procedure. To 
save computer memory in the larger lattices, we use a self-recurrent algorithm 
(in C language). In this way the total memory employed to sketch the clusters 
is almost negligible (it grows nearly as the lattice size squared). 

Due to the absence of MC dynamics, the system is specially vulnerable 
to eventual pathologies of the random number generator. We have observed 
significant deviations in some quantities for a commonly used shift register 
generator specially in the larger lattices. To avoid these effects, we have 



used as generator the sum (modulus 1) of the output of the generator of ref. 



TT| and a congruential one, since it is known that their respective drawbacks 
are very different.]^ 

To define the observables that we measure, it is useful to consider a related 
model that is a diluted Ising model with nearest neighbors infinite coupling, 
where the spins, (jj = ±1, live only in filled sites. It is easy to show that the 
magnetization of the latter model, 

i 

V being the volume, coincides with the strength of the percolating cluster in 
the thermodynamical limit and at T = 0. 

Knowing the size of the clusters, as their spins must take the same sign, 
we can write ^ 

c 

where Sc is the sign of the cluster c, ric its size, and the sum runs over all 
clusters. As Sc are statistically independent, we can construct an improved 
estimator for even powers of Ai (the only non-vanishing in a finite lattice) 
averaging over all possible values of {sc}, that henceforth we will denote as 

^ We have used Xn+i = 16807X„ mod(2'^^ — 1) for the congruential random generator, 
whereas the shift register formulas read: X„ = + using as pseudorandom 

number X„ XOR X„^6i- 
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(■ ■ •). For the second power we have 



■^^ = ^E^c- (4) 



c 



We define the susceptibihty as 

X = v(m^). (5) 

To compute the Binder parameter Vm we can construct an improved esti- 
mator for the fourth power of the magnetization. Averaging over signs, we 
obtain after some algebra 

c 

from which 

^ _ 3 1 {M^) 

For the finite-size scahng (FSS) method that we employ, it is very useful 
an accurate measure of the correlation length. We have used the second 
momentum definition [13| in the associated Ising model, that, in a finite 
lattice, reads 

X/F - ^ ^^/^ 
,4sin2(7r/L), 

where F is defined in terms of the Fourier transform of the magnetization 



r 



as 

F = - {||jW(27r/L,0,0,0)p + permutations^ . (10) 
It is also possible to construct an improved estimator for as 



^For another application of the Binder cumulant in percolation theory see jlj 
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To measure the critical exponents we use a form of the FSS ansatz that 
only involves measures on a finite lattice. For an operator O that diverges as 
{p — Pc)~^° 1 its mean value in a size L lattice can be written, in the critical 
region, as 

0(L,p) = L^o/^^ {Fom,p)/L) + 0{L-^)) , (12) 

where Fq is a scaling function and uj is the universal leading corrections-to- 
scaling exponent. From a Renormalization Group point of view, uj corre- 
sponds to the leading irrelevant operator. 

We can eliminate the unknown scaling function using the values from two 
different lattice sizes measuring at a p value where the ^/L quotients match. 
Specifically, defining 

Qo = 0{sL,p)/0{L,p), (13) 

we can write 

Qo\q^=, = s^o/^ + 0{L-). (14) 

Other examples of application of this method can be found in refs. [Q. 

The form of the scaling corrections allows to parameterize the finite-size 
effect on the determination of the critical exponents as 

^) -i"^) ocL-. (15) 

V / oo \ V / {L.sL) 

To compute the uo exponent, we can use equation (|T2|) for an operator 
with xo = (as, for instance, Vm or H L) obtaining for the shift of the 
crossing point of lattice sizes L and sL 



Ap^'^^ = [p,(L,sL) -p,(oo)] (X i^A^L— ^. (16) 

— 1 

To efficiently use the FSS formulas, it is necessary to use a reweighting 
method to move in the critical region. For this model there is not a Boltz- 
mann weight, but the role of the energy is carried out by the density of the 
configuration, and the probability distribution is binomial. 

The probability of finding a density q when filling sites with a probability 

p is 

From a set of measures of an observable O and the actual density of the 
configuration {(Oj, g^)} we can compute the mean value of the observable for 



6 



a neighbor density p' as 

Using equation (p!8D p-derivatives of observables can also be computed. 

Obviously we cannot extrapolate much further than \/p{l — p) /V, which 
is the dispersion of the distribution (0). Therefore the visible region de- 
creases as L~'^/^. Fortunately, it is enough for our purposes since to use 
eq. (0) we need to move in a neighborhood of the critical point whose size 
decreases as L~'^~^/''{^ L"^-^). 

3 Numerical Results 

We have produced a million of independent samples for each lattices, with 
L = 8, 12, 16, 24, 32 and 48. 

To measure the thermal critical exponent we have used as operators: 
dlogx/dp (xdiogx/dp = 1) and d^/dp (xd^/dp = I + u). For the magnetic 
exponents we have used the susceptibihty x i^x = l)- We remark that, 
although X is a fast varying function of p at the critical region (see refs. [§]), 
the use of eq. (|T^ allows a very precise measure. Moreover as what we 
directly measure is the quotient 7/z/ = 2 — 77, we can obtain a very accurate 
determination of the anomalous dimension 77. 

We have checked the method in the d = 2 case, where there is a very solid 
conjecture 0] for the values of the critical exponents, which is confirmed by 
conformal group analysis. We present the measured critical exponents for 
the two dimensional site percolation in table |l|, obtained from a million of 
samples for each lattice size. The conjectured values by Nienhuis |^ are 
7] = 5/24 = 0.20833 . . ., z/ = 4/3 and a; = 2. The agreement is very good. 

In the four dimensional case (see table |^), we observe a very stable value 
for the u exponent when using the operator d^/dp. However, the results for 
the exponents t] or u computed from measures of other operators do need an 
infinite volume extrapolation, what will be considered next. 

To measure the critical density and the corrections-to-scaling exponent 
u, we have studied the crossing points of Vm and ^/L for different pairs of 
lattice sizes, fitting the displacements to the functional form ([T6|) . As the 
behavior of Vm and ^/L is very different regarding the corrections-to-scaling, 
we obtain a great improvement performing a joint fit. 



7 



L 


V 


r] 


de/dp 


dlog(x)/dp 


X 


24 


1.324(9) 


1.326(14) 


0.2155(5) 


32 


1.330(8) 


1.30(2) 


0.2121(4) 


48 


1.344(10) 


1.36(2) 


0.2085(4) 


64 


1.330(9) 


1.36(2) 


0.2082(4) 



Table 1: Estimates for the critical exponents of two dimensional site percola- 
tion obtained from the finite-size scaling analysis using data from lattice sizes 
L and 2L. In the second row we show the operator used for each column. 
Practically we can read from the last row the conjectured values. 



We show in figure Q the crossing points of Vm and ^/L as a function of 
^ where we have used u = 0.689 and u; = 1.13. 

We fix the lattices ratio to s = 2 and perform the fit twice, for L > 8 and 
for L > 12. In both cases we obtain compatible values for the u exponent and 
for the critical density. We get acceptable fits, for example x^/d.o.f. = 4.7/4 
for the former. We give the central values from the former fit but with the 
error bars coming from the latter fit: 

u; = 1.13(10), pc(oo) = 0.196901(5). (19) 

The error bars have been slightly (20%) increased to take into account the 
error in the value of u. 
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Figure 1: Pc{L,2L) as a function of L ('^+1/'^) for the observables Vm and 



L 


V 


ri 


d^/dp 


dlog(x)/dp 


X 


8 


0.689(3) 


0.668(3) 


-0.0687(7) 


12 


0.687(3) 


0.666(4) 


-0.0775(7) 


16 


0.688(4) 


0.681(5) 


-0.0823(6) 


24 


0.691(5) 


0.683(6) 


-0.0868(8) 


00 


0.689(10) 


0.683(12) 


-0.0944(17+11) 



Table 2: Critical exponents obtained using data from lattice sizes L and 
2L for the four dimensional site percolation. In the second row we show 
the operator used for each column. The last row corresponds to the infinite 
volume extrapolation using (p!5|). 
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d 


UJ 


V 


7 


r] 


6 





0.5 


1 





5 


0.79 


0.57 


1.18 


-0.07 


4 


1.52 


0.68 


1.44 


-0.12 


3 


2.23 


0.83 


1.81 


-0.18 


2 


2.95 


1.07 


2.41 


-0.25 



Table 3: Our results using [2,1]-Pade resummation for the lo exponent. We 
also report (columns three to five) the results for the critical exponents iy^ 7 
and ri = 2 — 7///) obtained in ref. using the [2,l]-Pade-Borel resummation. 



Using these values, we can obtain an infinite volume extrapolation for 



the critical exponents by means of ([15|) . To control that higher order scaling- 



corrections can be neglected, we use an objective criterium. We perform 
the fit considering data from lattices of sizes L > Lmin and then repeat it 
discarding the smallest lattices data. If both fits parameters (extrapolated 
value and slope) are compatible, we keep the central values from the former 
fit and error bars from the latter. We have found that Lmin = 8 is enough 
for our data. 

The results are displayed in the last row of table 0. For t] the first term 
in the error have been obtained considering u fixed, and the second one 
corresponds to the variation when u moves within its error bars. In figure 
we show the behavior of t]{L,2L) as a function of L^'^, with u = 1.13, 
together with the extrapolated value. 

At this point we can use the results of ref. obtained with the e- 
expansion. We are specially interested in the corrections-to-scaling expo- 
nent, that is, the derivate of the /3-function at the non trivial fixed point (i.e. 
g* ^ 0). Using the /3-function and the non trivial fixed point from ref. we 
have obtained 



dm 

00 r- 



dg 



= e - 0.7607676^ + 2.00886e^ + O(e^). (20) 

g* 



As e(= 6 — d) is large, we have analyzed this series using the Pade technique. 
Only the [2,1]-Pade gives consistent results (i.e. u = 0(e) > 0). This agrees 
with the results of ref. where in the final Pade analysis of their series for 
the critical exponents only the [2,1]-Pade is reported (the results of the other 
Pades turned out to be incompatible with the numerical simulation results). 
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Figure 2: ri{L, 2L) as a function of L"'^, where lj is the corrections-to-scahng 
exponent, that we have fixed in this plot to 1.13 (see text for more details). 
The numerical value found in the literature, i] = —0.12, is also displayed. 

We show our results for u in table 3. We also display in this table the results 
for the exponents z/, 7 and i] calculated in ref. ||^ using the [2,l]-Pade-Borel 
resummation. 

In other cases with e = 2 a good agreement has been found between 
resumed series and numerical results. For the two dimensional Ising model 
the differences in rj and u are 1% and 5 — 35% respectively (taking into 
account the error bars of the e-expansion estimate of to')[0. However, we 
have obtained a discrepancy of 30% in the anomalous dimension and of 50% 
in the cj-exponent. Linking this discrepancy with the behavior of rj with 
the dimension, reported in the introduction, we find the e-expansion not 
trustworthy in this case. 
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4 Conclusions 



Using FSS techniques we have obtained accurate values for the critical expo- 
nents of the four dimensional site percolation. We have been able to param- 
eterize the leading corrections-to-scaling what allows to largely reduce the 
systematic errors coming from finite-size effects. 

We have obtained an anomalous dimension that is 30% far away from 
previous numerical and analytical (e-expansion) approaches. 

We project to extend these methods to the case of the diluted Ising model 
in four dimensions, in order to study the possible variation of the critical 
exponents on the critical line [^]. 
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